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We investigate a harmonically trapped two-component Bose-Einstein condensate within the miscible regime, 
close to its boundaries, for different ratios of effective intra- and inter-species interactions. We derive analytically 
a universal equation for the density around the different boundaries in one, two and three dimensions, for both 
the coexisting and spatially separated regimes. We also present a general procedure to solve the Thomas-Fermi 
approximation in all three spatial dimensionalities, reducing the complexity of the Thomas-Fermi problem for 
the spatially separated case in one and three dimensions to a single numerical inversion. Finally, we analytically 
determine the frontier between the two different regimes of the system. 


I. INTRODUCTION 

The experimental realisation of Bose-Einstein condensa¬ 
tion (BEC) in dilute atomic vapours was a landmark achieve¬ 
ment of late twentieth century physics [1-3], and there are 
now many elements that can be cooled down to quantum de¬ 
generacy [1-10]. A significant number of two-component 
Bose-Einstein condensates (TCBEC) have also been reported, 
as mixtures of two atomic species [1 1-14], two isotopes of the 
same species [15], or as two hyperfine states of the same iso¬ 
tope [16-23]. 

Many theoretical studies have addressed the density pro¬ 
files of TCBECs depending on the ratio between the intra- and 
inter-species interaction strengths [24-28]. However, most 
of these theoretical studies are numerical or based on the 
Thomas-Eermi (TE) approximation. This approximation, in¬ 
troduced for the single component case [29, 30], describes the 
basic features of the ground state of a BEC with large inter¬ 
atomic interactions. The TE approach neglects the kinetic en¬ 
ergy term in the time independent Gross-Pitaevskii equation, 
on the grounds that its contributions are to a significant extent 
dominated over by those due to the nonlinear interaction term. 
It can give good approximations, for instance, of the conden¬ 
sate chemical potential or of the order parameter near its max¬ 
imum value. However, close to the order parameter bound¬ 
aries, where the atomic densities are low, the TE approxima¬ 
tion cannot provide the condensate density profile. Knowing 
the wave function of the condensate around these boundaries 
is very important to characterize for instance the actual kinetic 
energy [31, 32], the tunneling rate across a potential barrier 
[31], or in the case of TCBEC systems, the penetration of one 
component into the other [27]. Several works have proposed 
new analytical approximations beyond the TE approach for 
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single component BECs [31-35] and for the two component 
case in the immiscible regime [27]. Here we present a new 
analytical approach to study the density profile of TCBECs 
within the miscible regime, around the regions where the TE 
approach fails, by deriving a universal equation. We also in¬ 
troduce a general procedure to solve the TE approximation 
of TCBECs in one (ID), two (2D) and three dimensions (3D) 
and we provide an analytical formula that determines the fron¬ 
tier between the different regimes of the system [24]. Our 
method also reduces the complexity of the numerical inver¬ 
sion required in the TE approach for the one and three dimen¬ 
sional cases [28, 36] for the spatially separated regime. 

The paper is organized as follows. In Sec. II we present the 
equations that describe TCBECs and the general form of the 
two-component TE approximation in the particular case of an 
isotropic harmonic potential. In Sec. Ill we derive a universal 
equation governing the behavior of the density profile close to 
the different boundaries of the system. Einally, in Sec. IV we 
develop a procedure to solve the two-component TE approach 
in a general way. We compare our approach, for each dimen¬ 
sionality, with the numerical solution of the coupled Gross- 
Pitaevskii equation in Sec. V, and present our conclusions in 
Sec. VI. 


II. GROUND-STATE OF A TWO-COMPONENT 
BOSE-EINSTEIN CONDENSATE 

A. Gross-Pitaevskii equations and the Thomas-Fermi limit 

The ground state of a TCBEC at zero temperature within 
the mean-field approximation is typically well described by 
the time-independent two-component Gross-Pitaevskii equa- 
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tions (TCGPEs): 

+ y, + gsN,\^'/ + -yu.j'I'. = 0 , 

( 1 ) 

where i = 1 or 2 refers to each component of the EEC, whilst 
iris, Ns, ViCr) and /ij are the mass, number of atoms, external 
potential and chemical potential of the s component, respec¬ 
tively. The intra- and inter-species interaction coefficients are 
given by gi > 0 and gi 2 , respectively. We will assume that the 
considered TCBEC is formed by atoms of the same species 
in two different spin states [16-23]. This means that we can 
set mi - m 2 - m. Eor simplicity, we consider equal trap¬ 
ping potentials yi(r) = y 2 (r) = y(r) which in our case will 
be isotropic and harmonic. Nevertheless, our results can be 
straightforwardly generalized for cases with mi + m 2 and for 
yi(r) y 2 (r). At this stage it is useful to redefine the interac¬ 
tion coefficients as gs - gsNs and gi 2 - gn V^ViAj, obtaining 
TCGPEs of the form: 

(-^ + V + 'I'., = 0 , (2) 

with n = A 1 /A 2 . 

By considering the TE limit, which neglects the kinetic en¬ 
ergy terms (V^T'j = 0) when compared with the nonlinear in¬ 
teraction terms, we are able to write down density profiles for 
either component. Adopting this limit the TCGPEs [Eq. (2)] 
become: 


2ms 


{v + gsfis + gnn^-^'Ws -iUs)'¥s^Q, (3) 

where we define iisir) - |T'j(r)|^ for s = 1 and 2. Then, by 
solving the two coupled equations (3) one obtains the general 
form of the TE density profile for each component in the re¬ 
gion where both components coexist, i.e. tis + 0, for both 
values of s\ 


(gi2n--3/2 _ g3_^)y(r) + 

glg 2 - g ?2 


(4) 

Note that in order to have positive-definite solutions within 
the TE approximation the denominator in Eq. (4), gig 2 - gj 2 , 
must be positive [37]. A system fulfilling this condition is 
commonly said to be in the miscible regime, otherwise it is in 
the immiscible regime. Throughout this paper, we will only 
consider intra- and inter-species interaction coefficients such 
that we are in the miscible regime. In the regions where one 
component is absent (n^ = 0 for i = 1 or 2) the density profile 
of the other component, within the TE approximation, reads: 


ns{r) 


P. - E(r) 

8s 


(5) 



R, Rs 


radius 

FIG. 1. Different density distributions obtained within the TF ap¬ 
proximation (solid lines) from Eqs. (4) and (5) and by numerically 
integrating (dotted lines) the TCGPEs [Eq. (2)] for the coexisting 
(a,b), and spatially separated (c) regimes. The different boundaries 
obtained in the TF approximation [Eqs. (6)] are highlighted. In this 
plot we assume gs > gis with If = 1. 


potential y(r) = V(r) = mw^r^/2, where air is the associ¬ 
ated angular frequency. With such a potential we can observe 
two different regimes: (i) the coexisting regime, where one 
of the components occurs only in coexistence with the other 
[Fig. l(a,b)]; and (ii) the spatially separated regime where 
both components occur partly in coexistence with each other 
and partly in isolation [Fig. 1(c)]. In general, we will denote 
the component with largest support (meaning the component 
with the largest spatial extent) with the subscript s. Within 
the TF approximation, we can distinguish two cases, when 
g 2 -s > gi 2 n''^^^ (equivalently g 3 _j > in which case both 
components have their density maxima at the center of the 
trap [Fig. 1(a)], or otherwise, when g 3 _j < g^TI*"^^^ in which 
case component 3 - s has its maximum of density at the cen¬ 
tre of the trap while component s has its maximum of density 
away from the centre [Fig. l(b,c)]. One can then note that if 
«s(0) > 0 (= 0) we are in the coexisting regime [Fig. l(a,b)] 
(spatially separated regime [Fig. 1(c)]). The condition sepa¬ 
rating these two regimes is derived in Sec. IV. 


B. Thomas-Fermi boundaries 

In a TCBEC, the external potential and the interaction pa¬ 
rameters determine the density distributions of the two com¬ 
ponents. In the following we consider an isotropic harmonic 


The Thomas-Fermi limit, as introduced above, can pro¬ 
vide us with a handle on the relative extent of each compo¬ 
nent. As such, we define the TF radii Rs, R^-s and Rs of the 
TCBEC trapped in an external harmonic potential by impos¬ 
ing iisiRs) = 0, iisiRi-s) = 0 in Eqs. (4) and iisiRs) = 0 in 
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Eq. (5), respectively: 

^ mu)l gl2n"-3/2 _ 
d2 ^ 2 gnUs - g 

mu1 gn - gin"-3/2 


(6a) 

(6b) 

(6c) 


In Fig. 1 we plot the density profiles of the TCBEC within 
the TF approximation and indicate the three boundaries given 
in Eq. (6); the outer boundary {Rs) that gives the maximum 
extent of the BEC, and the inner (Rj,^s) and innermost (Rs) 
boundaries delimiting the regions where the two components 
coexist. Note that the innermost boundary appears only when 
risiO) - 0 (i.e. R^ is only dehned when we are in the spatially 
separated regime). Figure 1 also shows the numerical solution 
of the TCGPE [Eq. (2)] in all the discussed cases, demon¬ 
strating that the TF approach cannot accurately describe the 
density prohles close to the boundaries. In the next section 
we will go beyond the TF approximation, deriving universal 
equations governing the prohles of the densities around the 
TF boundaries (the low density regions), thus enabling an im¬ 
provement on the approximate analytical prohles. 


III. UNIVERSAL EQUATION 

A. Overview 

In this section we present an analytical procedure to obtain 
the density prohle of a TCBEC around the outer, inner, and 
innermost boundaries. We derive a universal equation that 
describes the density prohle around each boundary, general¬ 
izing to two components the method developed in [31] for the 
single component case. This generalization will require addi¬ 
tional approximations for the inner and innermost boundaries, 
beyond those for the one component case. Note that in all the 
expressions shown in this section the subscript s refers to the 
component with largest support. 


B. Outer boundary 

In the vicinity of the outer boundary, Rs (present in both co¬ 
existing and spatially separated regimes), the TCBEC behaves 
effectively as if it had a single component. Therefore, we can 
follow the lines of [31] and linearize the harmonic potential 


around Rs'. 


V(r) - V(Rs) + mco%{r - Rs) + 0((r - Rs)^ (7) 


We then introduce this linearization into Eq. (2) with 'T 3 _ 
0, obtaining: 


fr j ~ 

^ + mojjRsir - Rs) + 


T', = 0 


( 8 ) 


where we have used that jj.s = V(Rs) [see Eq. (5)], and we 
have only kept the second derivative term of the radial part 
of the Laplacian, i.e., in two (three) dimensions we impose 
r^d'V/dr « « d^'¥/dr^). This approx¬ 

imation applies for values of Rs much larger than the thickness 
of the boundary given by Eq. (10), as discussed in [31]. 

Finally, by dehning the dimensionless variable 




r- Rs 


(9) 


with 


-L 


n 


,2 \l/3 


\2rrP-cjjRs 

and the dimensionless wave function cps through 


^s(r) = 


n 


ds ^J2mg, 




( 10 ) 


( 11 ) 


we obtain the following universal equation describing the pro¬ 
file of the outer boundary [31]: 




( 12 ) 


C. Inner and innermost boundaries 

Around the inner and innermost boundaries, Rs and R^-s, 
respectively [see Fig. 1(c)], both components coexist, and we 
must therefore consider the full coupled Eqs. (2). Thus, in 
order to obtain the density prohle around the s component 
boundary (the innermost boundary) we linearize the potential 
around Rs'. 

V{r) - V{Rs) + mc^Rsir - Rs) + 0{(r - Rs)^). (13) 

Introducing this linearization into Eqs. (2), one obtains the fol¬ 
lowing two coupled equations for the s and 3 - s component, 
respectively: 


2m dr^ 
2m dr^ 


- Rs) .,2„ I, ^ „ IVU ,2 , „..tt^-3/2,u,, ,2 


.gi2n 


s~i/2 


■83- 


+ mcjfRsir - Rs) + 


'P. =0, 


+ ma3%ir - Rs) + g3-.l'I'3-.l" + 


gn^ 


83-s 


(14a) 


'P 3 -. = 0 , 


(14b) 
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where the term V(Rs) in both equations has been rewritten us¬ 
ing the expression (6a). As for the outer boundary, we only 
keep the second derivative term of the Laplacian. The influ¬ 
ence of the first derivative is much less than that of the second 
derivative in the limit of large R^. Thus, this approximation 
will not be valid in the cases for which Rs is close to the ori¬ 
gin. Specifically, Rs must be larger than the thickness of the 
boundary [Eq. (17)] [31], or, in other words, the relationship 
between the nonlinear parameters has to be such that the sys¬ 
tem is far from the crossover condition between the coexisting 
and spatially separated regimes, for which Rs - 0 (analytical 
expression shown in Sec. IV). 

In order to solve the two coupled Eqs. (14) for the s com¬ 
ponent, we use the TE approximation for the 3 - s component 
by assuming that, close to Rs, the density of component 3 - s 
is large enough to ignore the kinetic energy terms, i.e., we 
impose d^'V^-sIdP' = 0 in Eq. (14b). There are two limit¬ 
ing cases where this assumption cannot be applied: (i) when 
R^_s - Rs - 0{e), which occurs when gi 2 —> yjgxgi, because 
the TE approach is at the limit of its applicability and (ii) when 
Rs - Ri-s = 0{e), which occurs for gi 2 —» 0, because there is 
no interaction between components and the system reduces to 
two noninteracting BECs. 

Then, by combining the TE form of Eq. (14b) with 
Eq. (14a), one obtains: 


2m dr^ 



^3-j 


ir-Rs) 


gs- 


g3-s 


I'P. 


T', = 0 


(15) 


By following a similar procedure one obtains the equiva¬ 
lent equation for the 3 - s component around R^-s at which 
n 3 _j(f? 3 _s) = 0. The resulting equation has the same form as 
Eq. (15) but with 3—s and s exchanged. Thus, in order to solve 
the innermost and inner boundaries we define the dimension¬ 
less variable: 


6 - ± 


r-Rs 

ds 


(16) 


with ds given by 


ds — 


1/3 


±2m^u)lRs (1 - 

and the dimensionless wave function (ps defined through 

n 




dsyj2m {gs-g]2/g3-s) 


U^s). 


(17) 


(18) 


In Eqs. (16) and (17) the + (-) sign applies for the inner (in¬ 
nermost) boundary, and for the inner boundary s has to be 
interchanged by 3 - i and 3 - shy s. 

Einally, one obtains the same universal equation derived in 
Sec. IIIB [Eq. (12)] for both boundaries: 


0';-(^, + 02)^, = o, 



- 4 - 3 - 2-101234 

FIG. 2. Numerical solution of the universal equation [Eq. (12) and 
(19)] (black solid line), and the two asymptotic behaviors given by 
the Hastings-McLeod solution [Eq. (20)]: the Airy function (green 
dotted dashed line), and (red dashed line). We also plot the 
asymptotic behavior of the Airy function [Eq. (21)] for ^ —» -l-oo 
(blue dotted line). 


where the s (3 - s) applies for the innermost (inner) boundary. 

Summarizing, in this section we have obtained the univer¬ 
sal Eqs. (12) and (19) that describe the density profiles at 
the boundaries of a TCBEC trapped in a harmonic potential. 
Note that we obtain the same universal equation for the three 
boundaries. However, they require different transformations, 
Eqs. (9) and ( 1 1) for the outer boundary and Eqs. (16) and (18) 
for the inner and innermost boundaries, in order to retrieve the 
actual wave function at each boundary. 


D. Solving the universal equation 


One can recognize Eqs. (12) and (19) in the literature 
[38, 39] as being a Painleve type-II equation, which for pos¬ 
itive defined solutions with no divergences or sinusoidal be¬ 
haviors has a Hastings-McLeod (HM) solution* [38] with the 
following asymptotics: 

( 20 ) 

1 for f ^ -oo, 

where Ai(^) is the Airy function. 

In this paper we are interested in the ^ » oo limit, where 

the density is small and the TE approach is not applicable. In 
this limit, the asymptotic formula of the Airy function gives 
an intuitive form of the behavior of the wave function at large 

^ +«,) ^ V2Ai(f) - (21) 


^ Note that the Painleve type-II equation has a factor two in front of the 
nonlinear term that makes the prefactor in front of the Hastings-McLeod 
solution equal to 1 instead of V2. 


(19) 
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Figure 2 shows the asymptotic forms of the HM solution 
[Eq. (20)] and the numerical solution of Eq. (12) or (19) for 
comparison. The asymptotic behavior of the Airy function 
[Eq. (21)] is also plotted. Even though the approximation 
given by Eq. (21) is obtained at ^ > oo, the two functions 

coincide even for low values of which makes this asymp¬ 
totic approximation very useful as an analytical expression to 
describe the universal equation for ^ > 0. It remains for us to 
find appropriate expressions for the TE boundaries in order to 
complete the transformations to retrieve the actual wave func¬ 
tions near each boundary. 


IV. THOMAS-FERMI SOLUTIONS 
A. Overview 

In this section we present a general procedure to obtain the 
radii and chemical potentials within the TE approximation of a 
TCBEC trapped in an isotropic harmonic potential in ID, 2D 
and 3D. We study both the coexisting and the spatially sep¬ 
arated regimes, and we determine the frontier between them 
for each dimensionality. 

The TE solution, either in the coexisting or in the spatially 
separated regime, can be found by following three general 
steps (a similar procedure was presented in [24]): (i) normal¬ 
ization of the density using the proper limits of integration, (ii) 
isolation of the chemical potential of each of the components, 
Us and as a function of the parameters of the system from 
the normalization integrals and (iii) insertion of /z, and 
into Eq. (4) [Eq. (5)] for the coexisting (non-coexisting) re¬ 
gion, to obtain the density profile. 

In some cases, these steps can be laborious, and in partic¬ 
ular step (iii) may not be available analytically. Specifically, 
we have found that in the ID and 3D cases the chemical 
potentials cannot be inverted analytically within the spatially 
separated regime. Thus, here we show how to reduce the two 
coupled algebraic equations to a single equation, reducing the 
complexity of the numerical problem [28, 36]. In all other 
cases, fully analytical expressions can be found. 


B. Coexisting regime 

The coexisting regime can be solved in a fully analytical 
fashion for all three dimensionalities. In this case, the nor¬ 
malization conditions for the s and 3 - v components read: 


«3-,, R, 

J" ({s + r?sr^)d°r + J' ( 


(^s + risr^) d^r + I (T, + x.r^) = 1, (22a) 

0 R 3 -S 

«3_-v 

(^3.s + T]3-sr^)d^r^l, (22b) 


/ 

0 


where D - 1, 2 or 3 depending on the dimensionality. The 
d^r differential represents the volume element for each case: 
d^r - 2dr, d^r - 2nrdr, d^r - Anp-dr, (note that in the ID 
case we add a factor 2 due to the fact that r should go from 
-00 to 00 ). We assume cylindrically and spherically isotropic 
configurations in the 2D and 3D harmonic potentials, respec¬ 
tively, and Tj, aTj, (s rjs are given by: 


As 


Ks 


gs^ 


m(x>l 
'2gs ’ 


glg2 - gi2 
gl2ff - g3_, 
2 glgl - gi2 


(23) 


By carrying out the integrations in Eq. (22a) and (22b) we 
reach: 


^oRi^s + eoRf = 1. 
Pd 

where Q^, and Po read: 


PdR^I^ = 1, 


Pd - -Aj 


2mul " - gs 


Qd — A; 


2mw; gi2 


- A/; 


3 gs 

2mcP: 1 


gig2 - gi2 
, glgl ~ gi2 3 


3 gs’ 

and where the scaling factors 
Ai = 1, A 2 = 


37T 

8 "’ 


A 3 - 


6ti 

15’ 


(24a) 

(24b) 


(25) 


(26) 


account for the different dimensionalities. Note that p^, is re¬ 
lated to D.O throughySo = -(gs/gi 2 U^~^^^)Q.o- 

Then, using Eqs. (24) and the definitions of Eqs. (6) we 
obtain: 






mo)^ 

^2/(D+2) 


n.-3/2^^ 


(1 - Q.DlPDf'^°^^^ , 

/ g, - _2/(0+2) 


(27a) 

(27b) 


Einally, by introducing these two chemical potentials into 
Eq. (4) and Eq. (5) one finds the solution of the TE density 
profile of a TCBEC in the coexisting regime. 

The analytical expression within the TE approximation of 
the frontier between the coexisting and spatially separated 
regimes can be found by using the fully analytical expres¬ 
sion of the density profile obtained by inserting Eq. (27a) and 
Eq. (27b) into Eq. (4) and setting ns(0) = 0: 


gl-s 


g\2 , 

-gl2(gl2 

gs 




-3/2 


gs) 


I gs + W-^'hn 
\g-‘^I^W-2l2g^2 



Note that Eq. (28) gives the condition that separates coexisting 
[Eig. l(a,b)] and spatially separated [Eig. 1(c)] regimes. All 
the calculations shown in this subsection are valid for D = 1, 
2 and 3. 
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C. Spatially separated regime 


As mentioned previously, the full solution of the two chem¬ 
ical potentials in the spatially separated regime cannot be 
found analytically in the one and three dimensional cases. 
Here we present a procedure to reduce the complexity of this 
numerical problem. We start by using the normalization con¬ 
ditions for the s and 3 - s components, respectively: 


Ri-.< 

„ 2 ^ jD 


I 

Rs 

Rs 

f 


+ 7jy)d^r+ J (As +^ 1, (29a) 

Ry-. 

R3-S 


(A3_s + K3^sr^)d^r+ J (^ 3 -s + r/ 3 -sr^)d‘^r ^ 1 , (29b) 

Rs 




,. 2 ^ jD,. 


where dj, Ks, ^s and tJs are defined in Eq. (23). After inte¬ 
grating Eqs. (29) and rearranging the terms we obtain the two 
coupled equations: 


+ ^dR^Is + ^dRT^ = 1 , (30a) 

+ PdR^s ^ I, (30b) 


Where 




ao - -A/) 


Imo)]. gi 2 

3 g^-s 


gl2 - g3-s 


,n3/2~ 


^1^2 - gn 


70 - Ad 


Imojj gnH" 2/2 


^ 1^2 




(31) 


with the scaling factors accounting for the different dimen¬ 
sionalities given in Eq. (26), and where Qd, €d and Pd are 
defined as in Eq. (25). 

In the ID and 3D cases the chemical potentials of both 
components cannot be obtained analytically from Eqs. (30). 
Therefore in order to reduce the two coupled equations into a 
single equation we rewrite and R^^j from Eq. (30) as: 

RT^ - ;- ^''^2 , [gsN' + , (32a) 

(^ 1^2 - gn^jD 

where A' = l-e^Rf^^- Introducing Eqs. (32) into the relation 
between the TE radii obtained from Eqs. ( 6 ) yields 


52 gsigZ-s - ^^^) o2 , ^12(^in" - ^12) o2 

K, - -^^-fi3_j. 


.gl^2 - ^12 


^l.g2 - gi2 


(33) 

and one obtains a single equation that only depends on jis and 
on the parameters of the system. This equation has two roots 
that need to be inverted numerically in order to find (the 
root power depends on the dimensionality). By using the ex¬ 
pression for /? 3 „j from Eqs. (32a) and ( 6 a) we find an analyti¬ 
cal formulation for p 3 _j as a function of /Zj. Once we have the 
two chemical potentials, we introduce them into the densities 


[Eqs. (4) and (5)], obtaining the TE solution of a TCBEC in D 
dimensions in the spatially separated regime. 

In the 2D case, however, Eqs. (30) can be solved analyti¬ 
cally using the definitions of a, (3, y, Q and e from Eq. (25), 
Eq. (26) and Eq. (31) for D - 2. In this case, using the ex¬ 
pressions of the TE radii from Eq. ( 6 ) we obtain: 


(2D) _ ^ / OT6j2g3^^(n3 2^ + 1) 

n 


„(2D) _ 


(34a) 


/ 23 -s + 


mw2(n" 2/2gj - gj2)(n" 2/2gi2-g3_,) 


7rn2'-3 


^12 


(34b) 


and by substituting into Eq. (4) and Eq. (5) we get the density 
profiles of the TCBEC in 2D in the spatially separated regime. 


V. DENSITY PROFILES AROUND THE BOUNDARIES 


In this section, we present the comparison between the 
density profiles in ID, 2D and 3D of a TCBEC around the 
boundaries defined in Eq. ( 6 ) obtained: (i) within the TE ap¬ 
proximation (Sec. IV), (ii) using the universal equation de¬ 
rived in Sec. Ill and (iii) by numerically integrating the TCG- 
PEs [Eq. (2)] for the spatially separated regime. To reduce 
the parameter phase space we consider the same number of 
atoms for both components, i.e., H = 1 , gi 2 = 3 lXgig 2 
gs - 2 g 3 -s, in such a way that the parameterdetermines the 
ratio between the different TE radii [Eq. ( 6 )], and g^-s deter¬ 
mines the strength of the nonlinear interactions. In addition, 
we rescale the densities of both components to the maximum 
value of the density of the 3 - s component at the origin (con¬ 
sidering the s component to be that with largest support) and 
the r coordinate to the maximum extension of the TCBEC 
(Rs)- We also use harmonic oscillator units, which is equiv¬ 
alent to setting H - m - air = 1. These settings allow us 
to compare the behavior, graphically, of the density profiles 
close to the boundaries for different values of the nonlinearity 
on the same axis scale. 

In order to have some reference values of the strength of 
the considered nonlinear interactions (g 3 _j) we compare our 
rescaled nonlinear parameters with typical experimental val¬ 
ues. We consider a TCBEC of trapped in a harmonic 
potential with a radial trapping frequency air -2nx 20 Hz, 
transverse trapping frequency aij_ = 27r x 150 Hz (for the one 
and two dimensional cases) and as = lOOao, with ao and as 
being the Bohr radius and the s-wave scattering length, re¬ 
spectively. In the ID case, g 3 _j ^ 200 corresponds to a BEC 
with Ns - N 3 -S - 1 X lO'* particles, while in the 2D and 3D 
cases, g 3 -s = 1000 corresponds to an approximate value of 
Ns - N 3 -S - 5 X 10"^ particles. 

Eigure 3 shows the density profile of a ID, 2D and 
3D TCBEC trapped in an isotropic harmonic potential in 
the spatially separated regime using the TE approximation 
[Eqs. (30)-(33)], the universal equation [Eqs. (12) and (19)] 
and the results of the TCGPE [Eq. (2)] for different nonlinear¬ 
ities. Note that we have fixed s = 1 for the component with 
















7 



FIG. 3. Comparison between the results of the density profiles obtained through the universal equation [Eqs. (12) and (19)], the TF approx¬ 
imation (Sec. IV) and the numerical solution of the TCGPEs [Eq. (2)]. The first row shows the TE density profile (blue dotted line) and the 
numerical solution of the TCGPEs (solid lines) of a TCBEC in the spatially separated regime for the ID (1st column), 2D (2nd column) and 
3D (3rd column) cases for two different values of g 2 in each plot. We also plot a magnification around Ri (2nd row), R 2 (3rd row) and (4th 
row) where we include the asymptotic approximation of the universal equation V2Ai(/-) (dotted-dashed line) derived in Sec. III. The parameter 
values used are: (i) for the ID case;^^ = 0.985, g 2 = 200 (black lines) and g 2 = 1000 (red/gray lines), (ii) for the 2D case;^' = 0.9, g 2 = 1000 
(black lines) and g 2 = 10000 (red/gray lines) and (iii) for the 3D case;^' = 0.9, §2 = 1000 (black lines) and g 2 = 10000 (red/gray lines). Note 
that the magnifications only include the component under study. Eor a complete description of the dimensionless parameters and scalings see 
the text. 
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largest support. We observe that in ID (1st column of Fig. 3) 
the asymptotic behavior of the universal equation at the dif¬ 
ferent boundaries is in excellent agreement with the numerical 
solution of the TCGPE, even for values of g 2 corresponding to 
a relatively small number of particles in a typical experimental 
TCBEC. 

In 2D (2nd column of Eig. 3) we can see that the universal 
equation around the boundaries (Sec. Ill) gives a good insight 
of the numerical solution of the TCGPE in a fully analytical 
way for g 2 - 1 x 10“*. However, for relatively small nonlineari¬ 
ties (g 2 ~ 1000), the universal equation close to the innermost 
boundary cannot describe the density of the TCBEC, as dis¬ 
cussed in Sec. 111. 

In 3D (3rd column of Fig. 3) we see that the outer and inner 
boundaries are in very good agreement for both nonlinearities, 
showing that this approximation can be used to describe the 
density of a TCBEC around the boundaries provided that the 
conditions mentioned in Sec. Ill are fulfilled. The universal 
equation close to the innermost boundary, on the other hand, 
hardly reproduces the density of the TCBEC for low values of 
g 2 , however, the approximation appears to have broad validity 
for values of g 2 above 10000. 

VI. CONCLUSIONS AND REMARKS 

In this paper we have presented an analytical approxima¬ 
tion to the ground state density profiles of a TCBEC trapped 
in an isotropic harmonic potential in the mean field approxi¬ 
mation around the boundaries of each component, where the 
TE approximation is no longer valid. We have derived univer¬ 
sal equations that give a very good estimation of the behavior 
of the density profile at the boundaries of each species, soft¬ 
ening the sharp edges produced by the TE approximation. We 
have compared our analytical results with the numerically in¬ 


tegrated TCGPE, obtaining an excellent agreement between 
them. The method proposed in this paper also offers the pos¬ 
sibility to calculate analytically, as proposed in [31, 32], the 
kinetic energy of the system, tunneling between double well 
potentials and other possibilities such as calculating an equiv¬ 
alent healing length in a TCBEC in the miscible regime, simi¬ 
larly to the penetration depth defined in the immiscible regime 
by [27]. Moreover, the approach presented in this work can be 
easily extended to different species (i.e. different masses), and 
also, due to its generality, the procedure may be extended to 
other trapping potentials. 

We have also studied the TE approximation for ID, 2D and 
3D. We have shown that the coexisting regime can be treated 
analytically in all three cases. However, the spatially sepa¬ 
rated regime only has analytical solution in 2D. In ID and 
3D we can decrease the complexity of the numerical inver¬ 
sion required by reducing the resulting system of two coupled 
equations to a single one. Einally, within the TF approxima¬ 
tion, we have determined, analytically, the frontier between 
the coexisting and spatially separated regimes. 
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